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We prove that when the equations are restricted to the principal part the standard version of the 
BSSN formulation of the Einstein equations is equivalent to the NOR formulation for any gauge, 
and that the KST formulation is equivalent to NOR for a variety of gauges. We review a family 
of elliptic gauge conditions and the implicit parabolic and hyperbolic drivers that can be derived 
from them, and show how to make them symmetry-seeking. We investigate the hyperbolicity of 
ADM, NOR and BSSN with implicit hyperbolic lapse and shift drivers. We show that BSSN with 
the coordinate drivers used in recent "moving puncture" binary black hole evolutions is ill-posed at 
large shifts, and suggest how to make it strongly hyperbolic for arbitrary shifts. For ADM, NOR 
and BSSN with elliptic and parabolic gauge conditions, which cannot be hyperbolic, we investigate 
a necessary condition for well-posedness of the initial- value problem. 



I. INTRODUCTION 

Numerical solutions of the Einstein equations are of- 
ten obtained by specifying initial data on a hypersurface 
and evolving in time. The Einstein equations then split 
into evolution equations which contain time derivatives 
and constraints which do not. A choice of variables and 
of evolution equations for them is called a formulation of 
the Einstein equations. At the same time, some metric 
components are not determined by the evolution equa- 
tions or constraints and must be determined otherwise to 
obtain a closed system. This is called a choice of gauge. 

A standard formulation is the ADM one, whose vari- 
ables are the metric 7^ and extrinsic curvature Kij of 
a spacelike hypersurface. The ADM evolution equations 
for these two tensors are first order in time and second or- 
der in space. The ADM constraint equations are second- 
order in space. The 7^- represent 6 out of 10 components 
of the metric g^^u of the spacetime. The other 4 compo- 
nents can be parameterised by a scalar a, the lapse, and 
a 3- vector /3% the shift. The lapse and shift do not ap- 
pear in the ADM constraints. They appear in the evolu- 
tion equations, but without time derivatives. They must 
therefore be determined by a gauge choice. 

One distinguishes fixed gauge choices, in which the 
lapse (or the densitised lapse) and shift are given a priori 
as functions of the coordinates x^ , and live gauge condi- 
tions. The latter can be subdivided into algebraic ones 
where the lapse and shift are algebraically related to the 
dynamical variables, differential (typically elliptic) con- 
straints on the lapse and shift, and evolution equations 
for them. 

In order for numerical solutions to converge to the con- 
tinuum solution, given that discretisation error in the 
evolution equations generates constraint violations, the 
continuum time evolution problem must be well-posed for 
arbitrary initial data that do not obey the constraints. In 



this paper we investigate the well-posedness of a number 
of live gauge conditions combined with the ADM formu- 
lation and three formulations that are derived from it, 
the NOR, BSSN and KST formulations. 



In Sect. ni we review the formulations while leaving the 
gauge still undetermined. We show that, roughly speak- 
ing, the NOR formulation is equivalent in the principal 
part to the BSSN formulation and the KST formulation. 
For the KST formulation this will require shedding three 
variables that evolve trivially. 



In Sect, mil we review live gauge conditions. Our guid- 
ing principle is that the gauge should be "symmetry- 
seeking" . Specifying the derivatives of the lapse and shift 
along the normal to the time slices makes it easier to es- 
tablish well-posedness analysis but specifying their time 
derivative makes it easier to show a gauge is symmetry- 
seeking. We show that the two time-derivatives are 
equivalent in some cases. 



In Sections llVIVII we combine formulations and gauges, 
and investigate the well-posedness of the resulting closed 
systems. What we can prove depends on the choice of 
gauge. In Sect . IIVI and Ivl we investigate gauges in which 
the whole system is strongly or symmetric hyperbolic. 
This gives us both necessary and sufficient conditions 
for well-posedness. In Sect. I VII we use mode analysis 
to find a necessary condition for the well-posedness of 
the Cauchy problem of the remaining, non-hyperbolic, 
systems. Sect. IVlTl contains our conclusions. 



II. EQUIVALENCE OF FORMULATIONS OF 
THE EINSTEIN EQUATIONS 

A. Definitions 

Well-posedness of an initial value problem implies that 
an estimate 



\\Sui;t)\\<FmSui-M\ 



(1) 



exists, where u(x,t) is the solution, u{x,0) the initial 
data, S denotes a linear perturbation and || • || stands 
for appropriate norms (which may involve spatial deriva- 
tives), and where F(t) is independent of the initial and 
boundary data. This means that the solution depends 
continuously on the initial and boundary data. Neces- 
sary or sufficient conditions for the well-posedness of the 
initial-value problem, such as strong hyperbolicity, can 
usually be determined from the principal part of the evo- 
lution equations. 

In this Section we discuss formulations of the 3-t-l Ein- 
stein equations derived from the ADM formulation with- 
out yet specifying the lapse and shift. These are systems 
of evolution equations which are first-order in time in 
all evolved variables, but second-order in space in some 
variables v and z and first-order in other variables w. 
We shall therefore be looking at quasilinear systems of 
equations whose principal part is 

(2) 
(3) 



V ~ dv + w + dz, 
w ~ ddv + dw + ddz 



where ~ indicates equality in the principal part, a dot 
denotes dt = 9/9t, and d stands for spatial derivatives. 
We will close this system later by providing evolution or 
constraint equations for the (gauge) variables z. 

We shall show that for certain parameter choices and in 
certain gauges these formulations are equivalent, roughly 
in the sense that there is a bijection between their vari- 
ables and their principal parts, and that the initial-value 
problem for one system is well-posed if and only if it 
is well-posed for the other. The detailed results will be 
stated in a Lemma at the end of each subsection. 



where Di is the 3-dimensional covariant derivative com- 
patible with 7ij, Cf3 is the Lie derivative along /?*, and 
Rij is the Ricci tensor of ^ij . Indices are moved implic- 
itly with 7y and 7*-' throughout this paper. (Note the 
distinction between =, which denotes the definition of a 
shorthand notation, and =, which denotes a constraint 
or, below, the definition of an auxiliary variable.) For 
convenience of notation we define the shorthands 



l^^lt],k, 



U 



7'Sjfe. 



(9) 



where a comma denotes a partial derivative. With these 
the principal part (second derivatives of 7ij ) of Rij is 



Ri 



'2 h^J-k' 



-ft, 



(^i) 



+ d 



iM)- 



(10) 



The ADM evolution equations are first-order in time in 
the variables v = jij and w — Kij. The highest spatial 
derivatives in the ADM constraints and the ADM evolu- 
tion equations are second spatial derivatives of jij and 
a, and first derivatives of Kij and /3*. The principal part 
of the Hamiltonian and momentum constraints are 



H 

M,, 






K_, 



(11) 
(12) 



Throughout this Article we consider the vacuum Ein- 
stein equations. The generalisation to matter is trivial 
in the typical case, for example perfect fluids or electro- 
magnetism, where the Einstein equations and the matter 
equations couple only through lower order terms. 



C. NOR 

The NOR formulation 2] is obtained from the ADM 
formulation by introducing the auxiliary variables 



/. 



P 



(13) 



where p is a constant parameter. This gives rise to the 
auxiliary constraints 



B. ADM 

We write the spacetime metric in the well-known 3+1 
form 

ds^ = -a^ dt^ + 7y {dx' + 13' dt) {dx^ + P^ dt) . (4) 

The ADM formulation of the Einstein equations in the 
form given by York U is 



da, J 

dtK,j 



Cp^ij - 2aK, 



ij, 



CfjKij - 

+a [Rij - ^-ivji^v J 
-2 



D.DjU 
2KaK' 



KK, 



H = R-K,jK'^ +K^ = 0, 
M, = DjK^,-D,K^O, 



(5) 

(6) 

(7) 
(8) 



G,; = .U-d, + |i, = 0. 

The full NOR evolution equations are defined as 
^a^J = ADM, 



dtK, 



't^zj 



ADM + aaG, 



dtfi = dt 



a, - -t, 



(«j')'' 



ADM 



-a7,j [cH + dGk' 
+ CpG, + 2baMi, 



(14) 

(15) 
(16) 
(17) 



where ADM in equations H15I16|) represents the right 
hand sides of IjSIBI) respectively, and the subscript ADM 
in (|17|l means that the time derivative of di and ti must be 
replaced by their expressions obtained from equation |(SJl. 
a, 6, c and d are constant parameters to be fixed later. 
(We have replaced the parameter a' of [3 by d = ca' .) 



In the following we use the shorthand 



do 



\dt-f3''dk), 



(18) 



to hide the transport terms in the principal part of the 
evolution equations. We also densitise the lapse as 



I7l 



t/2 



Q, 



(19) 



where Q rather than a is now considered as the dynam- 
ical variable. (Later we will often set tr = 0, so that 
Q = a. However, ct > is necessary to obtain a hy- 
perbolic formulation with fixed lapse and shift.) The 
principal part of the NOR evolution equations is then 



+ (1 - a)c^(jj) + -^{ap - 1 - o-)i(,j) 



(20) 



(21) 



(22) 



In contrast to the ADM evolution equations, the NOR 
evolution equations contain second spatial derivatives of 
the shift in the evolution equation for /*. We define v = 



7^J, 



(InQ,/?') and w = (K^jJ^). With a^b 



d — 0, the evolution of fi can be ignored, and ADM is 
recovered as a limiting case of NOR. 

The Hamiltonian constraint of the NOR system is de- 
fined only up to use of the constraint Gi. For definiteness, 
we denote by H the ADM Hamiltonian defined in (O in 
terms of jij and Kij but not ft. The evolution equa- 
tions imply a closed constraint evolution system whose 
principal part is 



doH 

doM, 



~2Mr\ 



(23) 



+ lG.y'-{l+2d)G,y^, (24) 



doG, ~ 2bM,. 



(25) 



In particular, the ADM constraint evolution is given by 
(I2:il24|l with a = b = c = d = 0. (The variable f, then 
decouples from 7,j and Kij and can be neglected). 



D. BSSN 

The BSSN formulation (see for example l4|) is obtained 
from the ADM formulation by introducing the new vari- 
ables 



% = (det7) ^/^7jj, 
I = 7-'7 7jfe,i, 



(26) 
(27) 



1 
12 



— Indet7, K = Y'K^j 



A, = (det7)-i/3 i^,, 



77u^ 



(28) 
(29) 



With a densitised lapse we then have In a = In Q + 6a(j). 
The definition of the F' gives rise to the differential con- 
straint 



G, = 7yF^ - 7^'=% 



?,Jfe^,. . , - 



0. 



(30) 



The definition of Aij gives rise to the algebraic constraint 



T = f^A, 



0, 



(31) 



and from the definition of 7^ we have the algebraic con- 
straint 



L»EElndet7 = 0. 



(32) 



(There are no definition constraints associated with the 
new variables (f> and K because "fij and Kij are no longer 
variables.) Here we do not define the full BSSN evolution 
equations but only their principal part. It is 



9o0 ^ 


6 6a 


(33) 


dojij ~ 


-2ly + 2a-' [p\,^lJ)k - ^%/3',fc) , (34) 


doK ~ 


-[(lnQ),,,+6a0,,,]e-4^r^' (35) 


doA,, ~ 


e-40 


-(lnO),», - 27""7»j-,mn - 2(1 + 3f7)0,,y 




+aik 


(^\o) + (1 - a)l'''lH^,l)l 


TF 

(36) 



9of' ~ 2{b-l)f't'A,k,i-ff'K,, 



(37) 



where TF indicates the tracefree part. (The entire right- 
hand side of (|36|l is made tracefree, not only the principal 
part.) The parameters a and b have been introduced for 
comparison with NOR. In standard BSSN, a = b — 1. 

For definiteness, we define the Hamiltonian and mo- 
mentum constraint for the BSSN system by substituting 
the inverse of the definitions H26|) and (|29|l into the ADM 
definitions of these constraints. We again use Mi and H 
to denote the resulting expressions. They do not contain 
the variables F\ Their principal part is 



H 

M, 



e-40^., {^fi^^^^^^ 



Arj^kl^'' 



:K_,,. 



(38) 
(39) 



The principal part of the evolution system for these con- 
straints is 



doH 



-2e-^'^f^M,, 



(40) 



9oM, c 


' ^^- + ^"''(1^''^-'= 






+ ^l''G,,>. + ll''D,,,), 


(41) 


doG, 2 


± 2&Mi, 


(42) 


doT c 


i 0, 


(43) 


doD c 


i -2T. 


(44) 



If T and Z? vanish initially, they vanish at all times, in- 
dependently of the other three constraints. (This is true 
for the full evolution equations, not shown here, as well 
as the principal terms given in H43I44|) .) From now on 
we restrict attention to the subspace of solutions where 
these algebraic constraints vanish. A discrete version of 
this restriction is implemented in the "BSSN" code of 
most groups by projecting to solutions of T = I? = at 
each time step. It is not clear how exactly this relates to 
the continuum restriction 5] . We shall call the restricted 
continuum system BSSN-C (C is for constrained). Note 
that although NOR has 15 variables and BSSN-C 17, two 
of those are algebraically redundant because we consider 
Aij to be tracefree and "fij to have unit determinant. 

There is a one-to-one correspondence between the vari- 
ables of BSSN-C and NOR given by (l2fil28l2Hll and 



E. KST with fixed gauge 

BSSN-C and NOR-B are both second-order formula- 
tions and have the same number of variables, namely 15. 
By contrast, the KST formulation ^ is first-order and is 
not evidently the reduction to first order of any second- 
order system J3|, while a first-order reduction of NOR 
has 33 variables, not 30. Nevertheless, in this Subsection 
we show that the two formulations are equivalent in the 
precise sense that when only the principal part is con- 
sidered, the KST formulation of the Einstein equations 
with densitised lapse and fixed shift is an autonomous 
subsystem, in 30 of the 33 variables, of a reversible re- 
duction to first order of NOR with densitised lapse and 
fixed shift. In the presence of an evolved lapse and shift, 
this equivalence must be revisited because the evolution 
equations for the lapse and shift could be written in a 
first-order form which is not a reversible reduction of 
their second-order form. This will be done in the fol- 
lowing two Subsections. 

In this Subsection, we shall assume that the lapse den- 
sitisation weight cr is a given parameter, and that the 
parameter p of NOR is also given. We shall prove that 
the remaining parameters (a, 6, c, d) of NOR correspond 
to the parameters (C, 7, x, rj) of KST as follows: 



From 



/,=7yr^' + (4-6p)0,, (45) 

we also have the useful relations 



^ l^''l^J,k 



40, 



U 



120, 



(46) 



The constraints H, Mi and Gi as defined for NOR in 
(I7l8ll4ll and for BSSN in i3xl3J)l3()ll are also equivalent 
in the principal part modulo T = D = 0, and we have 
anticipated this by using the same notation for both. Fi- 
nally, if c = --1/3, d = —a/3 and p = 2/3, the principal 
parts of the evolution equations and constraint evolution 
equations are also equivalent (for any a and b). The 
values of a and b in NOR then correspond to those of 
BSSN-C, and we have anticipated this by using the same 
notation for both. 

In the following we focus on two versions of NOR: one 
where a — b— l,c = d = Q, which wc shall call NOR- A, 
and one where a — b^l, c — d— —1/3, which we shall 
call NOR-B. ADM is also considered in the form of NOR 
with a = b = c = d = 0. In these three cases we also 
assume cr = unless specified otherwise. Our results can 
be summarised in the following 

Lemma 1: For any choice of lapse and shift, there is 
a one-to-one algebraic correspondence between the vari- 
ables of BSSN-C and of NOR-B with p = 2/3 and, using 
this, between the principal parts of their evolution equa- 
tions, constraints, and constraint evolutions. In partic- 
ular, BSSN-C is strongly/symmetric hyperbolic (in the 
definition of ^]) if and only if NOR-B is. 



46 = 2x + 477-p(3x + 77), (47) 

2c = 7, (48) 

4ab = -(1 + 3a)x + {l - '^ ^ ^c) V, (49) 

4db = -(1 + 27)x + (^ + ^C + 7) V- (50) 

Note that 5, x a-nd 77 only appear in the combinations 77/6 
and x/b- 

We define the shorthand 



ICij{Xklrn,C,'y,0-) 



1 + C 



X 



{i])k 



'IpiO Y^k 



2 

, 1 - c 

-\r k,l 

~ -^k ,1 



kij' 



X 



k {t,j) 



(51) 



for any Xkim that is symmetric in its last two indices. 
With the parameter identifications (|48I50|I , /Cy obeys the 
identity 

X V 

ICij{Xklm + -^JlmYk + -^lk{lYm)Xili<^) 

~ /Cy {Xkl„^, C, 2c, a) + aFf, j) + d-f.^Yk^'' (52) 

for any Xkim symmetric in its last two indices and 1^. 
Here ~ indicates that we have neglected first derivatives 
of jij. We now go from NOR to KST in three steps. 



1. First-order reduction of NOR 

Using our shorthand notation, we can write the NOR 
evolution equations with fixed gauge in the form 

9o7y - -2/^y, (53) 

^of^ ~ -2K,k'''+pK,, + 2bM,, (54) 

doKi-j ~ ICij (dk^im, *, 2c, a) + aGi^^j) + d-fijGk'%55) 

where the * indicates that the formal parameter C of /C 
cancels because jij ki — lij ik ■ Gi was defined above in 

We reduce the NOR system to first order by introduc- 
ing the auxiliary variables e^ij = lij,k- The new evolu- 
tion equations are 

di^Kij ~ JCij {ekhnXi 2c, a) + aGi^^j) + ^7^ G'fe''',(56) 
doCkij — ~2Kij^k, (57) 

with (|53|l and (|54(l as before. C now arises as a re- 
duction parameter, and we have used the shorthand 
Gi = fi — e'^ki — {p/2)eik'' for the new definition con- 
straint. Because no constraints have been introduced in 
(|57|l , the reduction is reversible in the sense defined in [3 . 
This means that the reduction is strongly (symmetric) 
hyperbolic if and only if the original second-order sys- 
tem is strongly (symmetric) hyperbolic. Furthermore, if 
both are symmetric hyperbolic, their energies are related 
by identifying ekij and 7y,fc. (See Q for definitions 
of strong and symmetric hyperbolicity of second-order 
systems, and the relation with their first-order counter- 
parts.) 

We now decompose Ckij (18 components) into ekij (15 
independent components) and Cj (3 components) as fol- 
lows: 



^kij 







Here H47|l is necessary for (|6()ll to hold. The two parts of 
Ckij evolve as 

Y Tl 

doCkij '^ doBkij - —jijdoek- —jkiidoCj), (61) 



X 
= efeij - ^lijBk- 


-^7/c(.e,), 


(58) 


k P k 
= e\,-'^e^k\ 




(59) 


- e^ ^p,^ 




(60) 



doCi 



-2K,j 



From (|5^ with Xkij - 
()56|l can be rewritten as 



^pK.^. 
Ckij and Yi 



(62) 
-e," we see that 



doKij ~ JCijiikhnX, 2c, a) + af^^i^^) + d7y/fc''', (63) 

The evolution equations of first-order NOR in final form 
are (|53I54I61I62I63() . When only the principal part is 
considered, the subsystem (15416 1I63() in the variables 
(s-kij, fi, Kij) is autonomous, and the variables (6^,7^^) 
follow passively, and so could, always in the principal 



part approximation, be found after the subsystem has 
been solved. Any estimate on the subsystem gives rises 
to a similar estimate for the full system. Therefore the 
full system is strongly hyperbolic if and only if the sub- 
system is. (We discuss symmetric hyperbolicity in Step 
3). 



2. From first-order NOR to KST 



We now define 



y Tl 

dk^j = ektj + ^ lijfk + T^ 7fc(Jj), (64) 



26 
X_ 
2b 



2b 

H 
2b 



SMj + 7^ lijGk + 7^ Ik(zGj). (65) 



Note dkij is equal to Ckij up to a constraint. From H61|) 
and (|54() , or equivalently from (|57|l and (|25|l , its evolution 
equation is 

dodktj ^ -2Kij,k + Xli3Mk + -n-ik(iMj). (66) 

e-kij and Yi -^ fi to 



Finally we use (|52|) with Xkij 
rewrite (|^ as 



df)K., 



'OJ^tj 



ICij{dklm,Cili'^)- 



(67) 



(|5^ . l|BB|) and ^7^ together form the KST evolution 
equations for fixed gauge 6j. (Note that a = 2i7kst and 
that in |g] the notation Dkij = dkij/2 is used.) In this 
second step we have discarded the dynamical variables e^ 
from the autonomous evolution system. Their place in 
dkij is taken by the fi. 

An equivalence between BSSN and a particular case 
of KST was established in Q without however clarifying 
the role of the algebraic constraints in BSSN, or noting 
the special role of the e^. 



3. Symmetric hyperbolicity 

The conserved energy of NOR with fixed gauge in its 
original second-order form is 

e{Kij,j^j^k,fi) = coeo + ciei +C2(d,^ -7.y,fc7*'''-') 

+C3[G, + 6(d,-t,)]', (68) 

where cqCq -\- ciCi is a complicated expression depending 
on (p, a, 6, c, d, a). (The special case c == d = is given in 
|3|]). The energy of NOR in first-order form is obtained 
by replacing 7^.^ with ekij, and setting C2 = (l-l-C)/4. In 
the first-order system C is a reduction parameter which 
then determines C2, while in the second-order system C2 
arises as a free parameter in the energy. 

The main difference between KST and first-order NOR 
is that the latter has the additional variable e^. However, 
there is a unique choice for C3 (a complicated expression 
involving (a, b, c, d, a) but not p) which, together with 
C2 = (1 -I- C)/4, eliminates e^ completely from the NOR 



energy. (This is non-trivial as three terms e*ei, e*/i and 
e^Cik^ must be eUminated from the energy with only the 
one parameter C3 to adjust.) The resulting NOR energy 
becomes the KST energy (compare App. A of ^) when 
Cfcij and fi are replaced by dj^ij and its partial traces 
according to (|64|l and with the parameter identifications 
and choice of C3 given here. Then KST with fixed gauge 
is symmetric hyperbolic if and only if NOR is. We 
summarize our results in this Subsection in 

Lemma 2: There is a reduction to first order (from 
now, NORl) of the second-order NOR evolution equa- 
tions (from now, N0R2) with densitised lapse and fixed 
shift. A subsystem of NORl (from now, NORla) is au- 
tonomous in the principal part approximation. It is possi- 
ble to relate the parameters of NORl and KST such that 
in the principal part approximation there is a bijection 
between the variables and evolution equations of NORla 
and KST. The evolution of the variables Ci which are in 
NORl but not in NORla is trivial in the principal part 
approximation. This means that NORl is strongly hy- 
perbolic if and only if KST is. This in turn means that 
N0R2 is strongly hyperbolic if and only if KST is. 

Furthermore, N0R2 admits an energy which, when 
translated into first-order form, does not contain Ci and 
so is an energy for NORla. There is bijection between 
this and the KST energy. This means that N0R2 is sym- 
metric hyperbolic if and only if KST is. 

(Here the hyperbolicity of second-order systems is the 
one defined in Q.) 



F. KST with evolved lapse: Sarbach-Tiglio 

The Sarbach-Tiglio (ST) formulation of the Einstein 
equations [13 is KST with the evolved lapse 



which obeys 



9o(lnQ;) ~ -fiLK, 



(69) 



where fiL = 2crcff — a~^dF/dK in the notation of [Tol| . 
We now show that the ST formulation is (an autonomous 
subsystem of) a reversible reduction to first order of NOR 
with this lapse. We start again from NOR and work 
towards ST. 

The live lapse NOR evolution equations are H53|) and 
(|54|l together with 

doK^j ~ -(lnQ:)^,y +/C(g.y,fc,*,2c,0) 

+aG(,j)+d7,jGfe'^ (70) 

We introduce the auxiliary variable 

Oi = (lna),j (71) 

and evolve it as 

d^Oi ~ -fiLK^i, (72) 

so that its introduction is reversible. In analogy with 
(|65|l we now define 



doA,~~fiLK^, + £.M,. 



(74) 



with ^ a new parameter. We complete the reduction to 
first order as in the fixed gauge case by introducing Ckij , 
then replacing its components e^ by fi to obtain dkij- 
This transforms (|70|l into 



doK^. 



'AiJ) +^('^fcjj,C:7,0)- 



(75) 



(ISSJ, (EEJ, (Ell and lEJ form the ST evolution equations, 
where the NOR and KST parameters are related by H47|l . 
(|15|l . (inOJl with (gnj replaced by 



Aab 



1 3 

2-2' 



CU-2t 



(76) 



Ai = at-\- —Gi 



(73) 



Compared to H49() . we have set ct = 0, and the addi- 
tional term —2^ appears on the right-hand side because of 
the appearance of (2f/6)G(i.j) in ^(ij). We have shown 
that ST is equivalent to a first-order reduction of NOR 
with the lapse evolution (|^ . The characteristic speeds 
squared of \\S\ simplify in our notation to Ai = /ii, 
A2 = 1 + 4(c — bd) and A3 — ah. 

Now hold \iL fixed and set (7 = 0. Then NOR has 
5 parameters (a, 6, c, d, p) and KST has 5 parameters 
(7i Cj Xi '7; C); but there are only 4 algebraic relations be- 
tween them, namely 2c = 7, and (for generic values of 
p) three linear relations between {b,ab,bd) and (x,?7,C)- 
One can consider p in NOR and C in KST as arbitrary in 
their system with no counterpart in the other evolution 
system. 

ST find a 3-parameter family of strongly hyperbolic 
formulations whose NOR equivalent is a5 = 1, d = ac. 
Their 2-parameter family is the special case c = —1/4. 
The 3 free parameters can be taken to be (7, ^, rj) in 
KST and {p, a, c) in NOR, with only 2 algebraic relations 
between these sets. ST also find a 2-parameter family of 
symmetric hyperbolic formulations, which corresponds to 
NOR with c = —1/3, d = —7/30. The free parameters 
can be taken to be (a, b, p) on the NOR side but on the 
KST side there are only two, for example (?7,C)i which 
are determined by (a, 6) alone. 

Note finally that the equivalence between NOR and 
KST would hold for any gauge that specifies the first 
time derivative of the lapse, with fixed shift. 

Lemma 3: Lemma 2 holds also for NOR and KST 
with a Bona-Masso type evolved lapse and fixed shift. 



G. KST with evolved lapse and shift: 
Lindblom-Scheel 



The formulation of Lindblom and Scheel (LS) 8] is 
again based on KST, but with an evolved shift as well as 
lapse. The gauge conditions are H111I112|) below, but the 
following comments apply to any gauge where the first 
time derivatives of the lapse and shift are specified. Be- 
sides the auxiliary variables d^ij and Ai already present 



in ST (called 2Dkij and Ti in LS), the auxiliary variable 
Mk = a~^l3^ ^f. is introduced, and all possible constraints 
with the right number of indices and level of derivatives 
are added to the right-hand side of its evolution equation. 
In our notation, these are 

d^Mk' = ... + i^iGk-' + V^5G\fc + ^64'Gr' + i>75k'H, 

(77) 
plus auxiliary constraints. (Our constants V'4...7 s-^e lin- 
ear combinations of V'4...7 of LS. The derivatives of Gi 
are added not in this explicit form but implicitly with 
the definition constraints of dkij and Ai. We think of Gi 
as the relevant part of those constraints, and the defini- 
tion constraints of ai and ekij as the part that is irrele- 
vant because its addition does not change the reversibility 
of the first-order reduction.) Because there are no con- 
straints whose time-derivatives are equal to these con- 
straint terms, we cannot use a trick analogous to H73I74|) 
to recreate the effect of these constraint additions in the 
first-order reduction within fully 2nd-order NOR. The 
smallest system fully equivalent to LS that could be con- 
structed would be a hybrid between second-order NOR 
and LS comprising the NOR variables {'jij, Kij, fi) and 
gauge variables {a, (3^) but also M'^. We have not anal- 
ysed this further. 

Conversely, the introduction of APk remains reversible 
if we set "04. ..7 = 0- In th^ notation of LS, this corre- 
sponds to the conditions 



16xCV'3 + (8 


x' + 8xv + W)i^4 


= 0, 


(78) 




■qipi - 4xV'5 


= 0, 


(79) 




■05 --06 


= 0, 


(80) 




^7 


= 


(81) 



on '04. ..7- With these restrictions, LS is a reversible first- 
order reduction of NOR with the same gauge. The iden- 
tification between NOR and LS parameters is then given 
by gTj), EHl), ISni) and IZH) and Vi = 2^. The remain- 
ing free parameters ('02, "03; "08: V'9: V'lo) of LS are reduc- 
tion parameters connected to the introduction of auxil- 
iary variables. They have no direct counterpart in the 
NOR evolution equations but can be identified with free 
parameters in the NOR energy. 

Lemma 4: Lemma 2 holds also for NOR with a 
Lindblom,-Scheel-type evolved lapse and shift and the sub- 
set of LS formulations characterised by the parameter re- 
strictions I78t81\}. 



III. ANALYSIS OF LIVE GAUGE CONDITIONS 

A. Symmetry-seeking coordinate choices 

Li the 3-f 1 approach the gauge freedom of general rel- 
ativity is fixed by specifying the lapse a and shift (3 in 
terms of the 3-metric 7^ and the extrinsic curvature Kij , 
and possibly given functions of (i,x*), through algebraic 
or differential equations. 



A unifying viewpoint for the classification of gauge con- 
ditions was given in [ij. A good gauge choice should 
have the property that if the spacetime is stationary the 
vector field [d/dtY = an"' -\- /3° becomes equal to the 
Killing vector, whatever the initial slice and the spatial 
coordinates on it are. We do not require that (9/9t)° be 
timelike but only that a = —na{d/dt)'^ > 0, where Ua 
is the unit normal on the surfaces of constant t. This 
is assumed in the following. In such coordinates, jij 
and Kij, as well as a and /3* become independent of t. 
Such conditions can be formulated in terms of the van- 
ishing of 3-1-1 time derivatives. Modifying the terminol- 
ogy of [ij we shall use the term "freezing conditions" for 
gauge conditions in which the metric is immediately time- 
independent in the presence of a Killing vector, and re- 
strict the term "symmetry-seeking" for gauge conditions 
in which the metric becomes time-independent asymptot- 
ically. 



B. Elliptic freezing conditions 

If the dynamical variables are frozen in the presence 
of a Killing vector, so are the lapse and shift. Freezing 
conditions for the lapse and shift can therefore not con- 
tain their time derivatives. An example is the iiT-freezing 
lapse 



K = 0. 



(82) 



For given /?*, jij and Kij, this is an elliptic equation for a, 
whose principal part is the Laplace operator associated 
with "fij . It is clear from the way we have written it that 
this is a freezing coordinate condition. Irrespective of the 
evolution system, we always use the simplest form 



K = P'K^, - Aa + aK.jK'^ 



(83) 



in any gauge condition. The principal part of this can be 
written as 



doK c:i -d,d'(\na). 



(84) 



li K = Q in the original slice (and therefore on all 
slices), the K-ireezmg lapse is called maximal slicing. 
Maximal slicing combined with minimal strain shift 



D'j^j = 0, 
or minimal distortion shift 

D%, = 0, 



(85) 



(86) 



were first suggested for use in numerical relativity in [13| . 

If, neglecting lower-order terms, we pull the time 

derivative out of the divergence, we can integrate these 

last two gauge conditions over t. In the context of the 

BSSN formulation, the gauge condition F* = has been 
considered under the name F- freezing shift. We consider 
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a generalisation in the form of the 2-parameter family of 
shift conditions 



where 



/. ^ 



di - -ti 



P = 0, 



sG, 



sfi + (1 - s)d^ H — U, 



(87) 

(88) 
(89) 
(90) 



where p and s are constant parameters. Essentially, /' 
is the divergence of the 3-metric. The price we pay for 
pulling the time derivative outside the divergence is that 
while jij transforms as a tensor under spatial coordi- 
nate transformations, /' does not. (However, p — and 
similar conditions (Dirac gauge) can be formally made 
covariant by writing /* as the covariant divergence of 7y 
with respect to a flat background metric Il2|). p parame- 
terises a conformal weight, and s — 1 gives us the option 
to express the gauge condition in terms of the NOR vari- 
able f, while with s = the gauge condition is expressed 
in terms of the 3-metric. With s = 1, p = p — 2/3 we 
have P = 7"i/3f i 

We shall loosely call all these gauge conditions F- 
freezing. They are elliptic in /3' for p < 2, and by con- 
struction are freezing. The principal part of /' = can 
be read off from 



doh ~ a-'j,k(3\/-' + (1 - p)«-i/3^',,, 
+2{sb - l)K,f' + {p-2sb)K^, 



(91) 



The principal part in /?* coincides with minimal strain 
shift and minimal distortion shift for p = and p ~ 2/3, 
respectively. 



C. Parabolic drivers 

Coordinate drivers turn elliptic equations for the lapse 
or shift into heat equations or wave equations by adding 
one or two time derivatives. Their motivation is that 
in a spacetime with a Killing vector, and with suitable 
boundary conditions, we expect that the solution of the 
resulting heat or wave equation will tend towards a (time- 
independent) solution of the original elliptic freezing con- 
dition. 

Parabolic drivers for the lapse and shift were, to our 
knowledge, first suggested in [lj|. The basic (explicit) 
parabolic "ii'-driver" is 



(In a)' — —vlK 



(92) 



for some constant V]^ > 0. The choice of powers of a here 
and in the following is determined by the requirement 
that the gauge choice be invariant under the rescaling 



of i by a constant factor. The equation is also invariant 
under arbitrary changes of the spatial coordinates. 

The elliptic equation K — Q has been turned into a 
heat equation for a. If all variables other than a were 
time- independent, the solution of the parabolic equation 
would tend to a solution of the elliptic equation K = 
0, but in reality all variables evolve together. Whether 
the solution actually tends towards a freezing condition 
will therefore depend on the physical problem. However, 
we can see that the driver is at least compatible with 
being symmetry-seeking in the sense that the evolution 
of a stops if and only if the freezing condition K = is 
obeyed. 

Integrating explicitly over t, we obtain the basic im- 
plicit parabolic _?C-driver 



Ina = —vlK + 4>{x). 



(93) 



As -ftT is a dynamical variable, this is now an algebraic 
gauge condition rather than a heat equation, and so the 
implicit and explicit forms are not equivalent as PDE 
systems. However, if the constraints are obeyed the two 
forms will generate the same gauge. In particular, 4>{x) 
allows us to still set gauge initial data for a independently 
of the value of K given by the geometric initial data, as 



(j){x) = {vlK + \na\ 



(94) 



(In the remainder of this Section the suffix denotes 
the initial value of a quantity.) On the other hand, if 
the parabolic driver succeeds in driving K to a time- 
independent value, 4>{x) also characterises a condition on 
the asymptotic values of K and a, namely 



[x) — {vlK + \\ia)c 



(95) 



where the suffix cxd denotes the asymptotic value. Such 
a condition is to be expected: while elliptic conditions 
freeze K and /* at whatever value they have in the initial 
data, a driver condition freezes K and /' only asymptoti- 
cally and at values which depend not only on the geomet- 
ric initial data Kij and 7^ , but also on the initial values 
of a and /?* (and a and /?* in the case of a hyperbolic 
driver, see below). 

The situation for the F-freezing shift conditions is very 
similar to the /f-freezing lapse condition. The basic 
parabolic "F-driver" is 



/3' = ^5(«f)- 



(96) 



for some constant vs- The left-hand side is a vector and 
the right-hand side is not (which means that the time evo- 
lution does not commute with generic spatial coordinate 
transformations), but the equation transforms correctly 
under rescaling the spatial coordinates by a constant fac- 
tor and under rotations. The factor of a makes the equa- 
tion invariant under rescaling t by a constant factor. For 
p = this is a heat equation for /?% and for p < 2 it is 
parabolic. Again, this condition is compatible with being 



symmetry-seeking. Similarly, an implicit parabolic shift 
condition is 



P' = lysaf-fix). 



(97) 



What would happen if we replaced dt with do? The 
straightforward replacement, 

doi\na) = -:yLdoK, (98) 

can be integrated using the method of lines: 



where 



Ina + vlK — (t)[x{x, i)]. 



x^=(3^x'„ x'{x,Q) = x' 



(99) 



(100) 



and (j){x) is given by (|Mjl . Clearly this couples the lapse 
and shift conditions through i*, leading to a gauge con- 
dition which is different from the original one. However, 
if we write 



doilna) = -vlOoK + a'^ f3' d^(j){x) , 



(101) 



this can be rewritten as the do derivative of H93(l . If 
we again define <j){x) from 1)94(1 . we have (|93|l identically. 
Taking a dt derivative, we obtain (|92|l . Therefore H92() 
and Hl()l|l give rise to the same gauge (when all Ein- 
stein equations are obeyed). In particular we have l|^ . 
However, the do and dt implementations have different 
principal parts, and so their well-posedness is not obvi- 
ously equivalent. The use of do in the gauge conditions 
simplifies the well-posedness analysis, but do (3^ is simi- 
lar to the nonlinear term in Burger's equation, and may 
develop shocks [l^- Contrary to the claim in Q, this 
problem cannot be fixed simply by introducing an auxil- 
iary variable for the first derivatives of the shift, as the 
reduction admits the same blow-up solutions. 



D. Hyperbolic drivers 

Hyperbolic drivers were discussed in 16]. The basic 
explicit hyperbolic ii'-driver is 



(a"i(lnay)'+2KL(lnaJ: 



'^J'L 



K. 



(102) 



for some /i^ > and hl > 0. This equation has the same 
transformation properties as the parabolic lapse driver. 
The motivation is that the elliptic equation K — will be 
solved by turning it into a wave equation with a friction 
term for a. Once again, this would be true if all variables 
other than a were time-independent but it is not clear if 
this gauge condition is symmetry-seeking in general when 
it is applied to a situation where all dynamical variables 
evolve. 

As we have chosen the outermost derivative to be dt, 
we can again integrate in t to obtain the basic implicit 
hyperbolic iiT-driver 



which is again equivalent as a gauge condition to the 
explicit driver as long as the constraints are obeyed. 4>{x) 
can be used to set initial data for a independently of 
a and K, but if the spacetime admits a Killing vector 
and the driver is indeed symmetry-seeking then it is also 
related to the asymptotic values of a and K: 



(f){x) — {fiLK + 2kl In a)c 



(104) 



Another damping mechanism (suggested by Shibata 
|17| for the shift driver) for an implicit hyperbolic driver 



IS 



-/i^a K — 2i'LaK 



a 



(105) 



By linearising in a, holding all other variables fixed and 
carrying out a Fourier analysis, one sees that all modes 
are damped. However, this linear equation for a (with 
all other variables held fixed) is neither hyperbolic nor 
parabolic. Shibata chooses 77 ex Ax, so that in the con- 
tinuum limit the implicit hyperbolic driver is recovered. 
One advantage of this damping mechanism is that now 
(j){x) — ^llKoo, giving better control of the final gauge. 

Again, a similar analysis applies to the F-freczing shift 
conditions. The basic explicit hyperbolic F-driver (with 
damping) is 



(a-i/3')-+2Ats/3' = (/i5af)-. 



(106) 



(This equation has the same transformation properties 
as the parabolic shift driver.) For p — Q this is a wave 
equation for (3^. The equivalent implicit hyperbolic F- 
driver is 

/j* = fisa^r-2Ksa(3'-a<j>'{x), (107) 

and the alternative damping mechanism is 



p' = ^^sa^^ + 2lysaf-a<j>\x), 



(108) 



the latter obeying ^'(a;) = ^sctflo- 

Finally we note that the explicit parabolic or hyper- 
bolic drivers differ from the explicit ones in that in the 
expressions for K or /* that appear on their right-hand 
side one can add suitable constraints in order to change 
the principal part of the equations, while this is not pos- 
sible in the implicit case. 

We can strip the outer time derivative off the ex- 
plicit hyperbolic drivers to obtain the implicit hyperbolic 
drivers even if that outer time derivative is do with the 
same trick of adding a source term we used above in 
(|101|l . This affects the nature and well-posedness of the 
system, but leaves the gauge choice unaffected. This is 
not true if we replace the remaining (inner) time deriva- 
tive dt with i9o on the left-hand side of the implicit hy- 
perbolic drivers. For example, with 

(9o(lna) = -/iLi4:-2KLlna + 0(x), (109) 

(|104|l is replaced by 

(j) = {^1LK + 2KL\na + a-"^ (i'a.^)oo, (HO) 



a = —a fiLK — 2KLa In a -I- a (j>{x), (103) which couples the asymptotic lapse and shift. 
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E. Boundary conditions 

Elliptic gauge conditions guarantee that 4 dynamical 
variables are frozen in any spacetime, but 7^, i^y-, a 
and /3* are all frozen in a stationary spacetime only if the 
boundary conditions are compatible with the Killing vec- 
tor. Finding boundary conditions which are at the same 
time compatible with the symmetry-seeking gauge, the 
constraints, and give rise to a well-posed initial-boundary 
value problem is an important problem that has not yet 
been studied in any depth. The parabolic and hyperbolic 
drivers require the same type of boundary conditions as 
the elliptic freezing conditions, namely Dirichlet, Neu- 
mann or mixed boundary conditions for the lapse and 
shift. These cannot be completely homogeneous, as they 
must fix an overall common factor in the lapse and shift 
(corresponding to the freedom to rescale t). We shall not 
otherwise investigate this issue here. 

Finally, we note that when the the gauge evolution uses 
do and a symmetric hyperbolic evolution system is ob- 
tained, one can arrange the characteristic speeds so that 
singularity excision surfaces inside black holes become 
purely outflow, and no boundary conditions are required. 
When the hyperbolic gauge evolution uses dt , or when the 
gauge condition is elliptic or parabolic, boundary condi- 
tions for the gauge are required at excision boundaries. 
This may be desirable or undesirable. From the fact that 
the Einstein equations are well-posed in harmonic gauge, 
where excision surfaces are purely outflow, it is clear that 
no physical information can leave the black hole even if 
boundary conditions are imposed inside. 



IV. HYPERBOLICITY WITH IMPLICIT 
HYPERBOLIC LAPSE AND SHIFT DRIVERS 

A. Generalised Lindblom-Scheel gauge 

In the following we investigate the strong and symmet- 
ric hyperbolicity of ADM, NOR-A and NOR-B coupled 
to a family of implicit hyperbolic driver gauge conditions. 
These conditions are 



{do + Qsdp)^' 



-fiLK + eLa-^(3\^, (111) 
fisaf + esa{\nay\ (112) 



where we have not written the friction damping terms 
and the integration constants 0, 0' because they are non- 
principal. We have introduced the shorthands /3* = 
a~^/3* and 9^ = (3^di. 

The free parameters in the principal part of the gauge 
conditions are p, (l, Cs, fJ-L, ^J■Sl ^l and eg. With 
Cl = Cs — these are the conditions considered by Lind- 
blom and Scheel [3 and the parameters /il, l^s, ^l and 
eg have the same definition. The parameter p is related 
to the parameter A of [g by p = 1 -t- A. Cl = Cs = sim- 
plifies the analysis of hyperbolicity, while Ci = Cs = 1 
corresponds to the hyperbolic drivers we have discussed 



above and gives rise to simpler asymptotic gauge condi- 
tions. Harmonic gauge is the special case Ci = Cs = 0; 
PL = f^s — 1, ^L = 0, es = —1, p = 1- We shall refer to 
(|111I112() in generality as LS gauge. 

The parameter s appears implicitly in the gauge con- 
ditions, but its value does not affect the gauge as long 
as the constraints are obeyed, and it should therefore 
be considered as a parameter of the formulation. Both 
NOR-A and NOR-B have the free parameter p, but in 
the following examples neither the characteristic speeds 
nor the well-posedness depend on the value of p. 



Using do 



4 ^ 


E A±B, 


(113) 


2A E 


E ^i + (2 - p)ps + ^se-L-, 


(114) 


A^ E 


E {2- p){€L - PL)fJ.S- 


(115) 



We now follow |2j and investigate (|111I112|I with ^^ — 
(s = 0, thus using do as the time derivative. In the princi- 
pal part, frozen coefficient approximation the equations 
decouple into three blocks, which we consider in turn. 
The statements below hold for ADM, NOR-A and NOR- 
B, unless further qualified. 

The transverse traceless symmetric tensor block of the 
principal part is always diagonalisable, with A = ±1 
(twice) . The transverse vector block has eigenvalues A — 
±1 (for NOR) or (twice, for ADM), plus A = ±y^. 
It is diagonalisable for ps > in the ADM case, and for 
ps > with either ps y^ 1 or s — 1 in the NOR-A and 
NOR-B cases. 

The scalar block always has A — 0, ±1, ±u±, where 



B' 



v± depend only on the gauge choice, not the parameters 
p, s of the formulation. Defining B to be positive, they 
are real if and only ii < B < A. At generic points in the 
interior of this region the scalar sector is diagonalisable 
because the four gauge speeds are different. On the lines 
A = B, B = and A - B = 1, A + B ^ 1 some of the 
characteristic speeds are repeated, and generically these 
marginal cases are not diagonalisable. In particular, A — 
B (i.e. w„ = 0) happens for p = 2 or ps = or cl = Pl 
and then the system is never diagonalisable. There are 
further parameter choices where some of the marginal 
cases become diagonalisable, but we do not analyse them 
here. We just mention that under the asumption A j^ B 
and using the fact that for ADM, NOR-A and NOR-B we 
have c — bd= 0, diagonalisability of the whole 7x7 scalar 
sector with eigenvalues (0, ±1, ±v-,±v+) is equivalent to 
diagonalisability of the matrix 

/ eLes + {2-p)ps -el 0\ 

-pLes-{2-p)ps PL (116) 

V -PL^S - (26s - p)ps 1 - 2a6 + pL 1 / 

with eigenvalues {l,vt,v'^). 

The results simplify if we partially decouple the lapse 
and shift conditions by setting e^ = 0. Then v± are 
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simply ,/]xh and \/(2 — p)/^s and the matrix (|116|) be- 
comes triangular, which allows a simple discussion of its 
diagonalisability: as we said, we always need /xl > 0, 
/is > 0, p < 2 to avoid a zero eigenvalue w^. If the three 
eigenvalues are different then diagonalisability is guaran- 
teed for both ADM and NOR. If /^l = 1 then we need 
ah = 1, which is only possible for NOR; if (2 — p)\xs = [i-L 
we need es = —1; finally if [2 ~ p)iis — 1 we need 
ifJ-L^s + 1)(1 — a6) = fis{bs — l)(/iL — 1), which requires 
bs — 1 for NOR-A/B. It is possible to make all scalar 
speeds equal to ±1 using the intersection of those three 
cases: eL — 0, es — —I, ah = bs = ^l — 1, c — bd and 
(2 — p)fis — Ij while s, d, p, p, a still are free parameters. 
We can set all speeds (including vector speeds) equal to 
zero or one by further setting 11$ — s — I, which then re- 
quires p = 1. The result is the implicit form of harmonic 
gauge in NOR with a = b = 1 and c ~ d, while d, p, a 
are still free parameters. For example, setting d = —1/3, 
p — 2/3 and cr = we obtain a system equivalent to 
BSSN-C with implicit harmonic gauge. (Note that al- 
though the gauge is harmonic gauge, this system is not 
equivalent to the usual harmonic evolution system.) 

In all other cases we cannot make all speeds equal to 
zero or one. This is not necessarily bad: one may want 
constraint- violating modes to have speeds less than 1 so 
that they do not pile up on the black-hole horizon [a| , and 
one may want gauge modes to have speeds larger than 
1 so that one can impose boundary conditions on them 
at the black hole excision boundary (inside the horizon) 
in order to control the gauge better [l^. We have not 
investigated these issues here. A simple example is 



p = 2/3, PL = PS 



1 

4' 



^L^^s-- 0, (117) 



which has speeds A = (0, ±1, ±1/2, ±1/V3) < 1, inde- 
pendently of p and s. 

Finally we note that the lapse driver alone, with fixed 
shift, is strongly hyperbolic with NOR-A/B for p^ > 0, 
but not with ADM. 



The fact that this condition depends on /3„ = /3^uji/\uj\ 
raises a serious problem: As long as \P\ = \/WPi > 
I /is — 1 1, this inequality will be violated for all uji that lie 
on a cone in wave number space. In particular, it is then 
violated for arbitrarily large \uj\. This means that for 
sufficiently large shift, the system is only weakly hyper- 
bolic. The inequality required for strong hyperbolicity is 
therefore 



<Im5-i| 



(119) 



Numerical evolutions should be unstable when it is vio- 
lated. However for s — 1 pi9|l is not required. 

The characteristic speeds in the scalar sector are A = 
0, ±1 and the four roots of the quartic equation 



A^ + {CL+Cs)PnX^ 

+ [ClCsPI - eLes + {p- 2)/is 

+ [(p-2)a/iS-CsML]^nA 

+ (eL-/iL)Ais(/O-2) = 0, 



Mi 



(120) 



For C-L = Cs = this is a biquadratic equation and its 
solution is given in (|113llll5|l . For (^ — (g = 1 the roots 
are more complicated. If we restrict to e^ = 0, we find 
that they are 




+ (2-p)/is. (121) 



Diagonalisability of the scalar sector of ADM and NOR- 
A/B with LS gauge is guaranteed in the generic case 
where (with e^ — 0) 



C. Using dt 

We now set Cl = Cs = 1 1 which is equivalent to replac- 
ing do by a^^dt as the time derivative on the left-hand 
side of I|111I112|) . As before all speeds A are with respect 
to normal observers and in units of the speed of light. 

The characteristic speeds in the vector sector are A = 
±1 (for NOR) or (twice, for ADM), plus 



(118) 



With ADM, the vector sector is diagonahsable for ps > 
0, for any /3„. With NOR, ps > is also necessary, but 
not sufficient: two speeds which are generically distinct 
coincide for /is — 1 = ±/3n, and we need to impose ps — 
1 y^ ±/3n to maintain diagonalisability. 




(2- 


P)PS - PL 7^ 0, 


(122) 




Iml-1| > l/3|, 


(123) 


1(2- 


-p)ps-l\ > \(3\ 


(124) 



(Note that the equivalent conditions for the do case are 
obtained by setting /?* to zero here.) However, with 
ab = 1 condition (|123|l is not necessary, and (|124|l is 
not necessary for NOR with s = 1. 

In conclusion NOR-A/B with LS gauge can be made 
strongly hyperbolic with s = 1 for all values of /3*, at 
least for e^ = 0. 

Finally we note that the lapse driver alone, with fixed 
shift, is strongly hyperbolic for NOR-B, but not for ADM 
(the vector sector fails to be diagonahsable) or NOR-A 
(the scalar sector fails to be diagonahsable. The same re- 
sult was found in the ST system [10|, which is equivalent 
to NOR with the lapse driver and fixed shift. 
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D. Symmetric hyperbolicity 

1. NOR with live lapse and fixed shift 
The most general conserved energy has the form 

5 



E' 



with 



£0 = K,^K'^ 
ab~ 1 



l^J.k' + d'Y, 



■ dC-d, 



hd 



■ ft, 



(125) 



{ad'+dt')X, (126) 



ei 



if^ + f [Y^ - (a + 3d)X, 
1 



-(2 - 2a6 + 6c + &bd + /xl)*, 



(127) 



£2 


= l''{dd,-l'''W'^l^k,^Jm,n). 


(128) 


£3 


- x,x\ 


(129) 


£4 


= x«r\ 


(130) 


£5 
ier( 


= Y,Y\ 
3 we have defined the shorthands 


(131) 




x^ = u~{b~^% + {b-i)d,, 


(132) 




Y, = (logQ),, + -(ff-/Xi)i,. 


(133) 



(This energy is found by writing down the most general e 
and flux </)* quadratic in 7y,fc, Q.i, Kij and ji constructed 
using only these and 7*^, and restricting their 14+10 free 
coefficients by 19 conditions arising from energy conser- 
vation e = 0' j in the high-frequency approximation. We 
do not give the fluxes here. Positive definiteness of e 
has not been imposed yet.) Energies ea, £4 and £5 are 
automatically conserved with zero flux. Conservation of 
Co£o + ciei requires 

co(l-2a6 + ^L + 2c-26d)-|-2ci(l-a6 + 3c-3M) = 0. 

(134) 
with non-zero flux, and conservation of £2 (with non-zero 
flux as well) poses no new condition on the parameters. 
The coefficient cq must be strictly positive and therefore 
there are two possibilities: either ab — \ ^ 3(c — bd) 
and then c\ is determined by cq and the parameters 
{a,b,c,d,^j,L), or we have the special case ab — 1 — 
3(c — bd) = 3{fiL ~ l)/4 and then cq and ci are inde- 
pendent. Both NOR-A/B belong to that special case 
and require /x^ = 1, which is harmonic slicing. Generic 
NOR however allows for arbitrary values of fi^ ■ 

With the evolved lapse we obtain the same slicing as 
with fixed densitised lapse if we set /^l = ct and dtQ = 0. 
Then Yi can be dropped from the energy and H125H134'|) 
reduces to the energy of NOR with densitised lapse and 
fixed shift (EHIl- 



Fixing an overall factor by setting cq ~ 1, the energy 
has cither 9 free objects (4 coefficients Ci plus 5 parame- 
ters), or 8 free objects in the special case (5 coefficients 
plus 3 parameters). The ranges for those objects are re- 
stricted by the positivity conditions on the energy. The 
tensor and scalar sectors require 



1 
ci>-3, 



1 



1 



-4<C2<2, 



(135) 



and the vector sector requires positivity of a 4x4 matrix, 
which is difficult to convert into explicit inequalities for 
the coefficients and parameters. Some necessary condi- 
tions are 



C3 >0, 



C5 >0, 



C4 < 4C3C5 



(136) 



For example for NOR-A with /j.^ = 1 and choosing C4 = 
we still have 4 free coefficients, and a complete set of 
symmetric hyperbolicity conditions are 



Ci > 0, < C2 < i C3 > ci -t- ji 



4C2' 



^ (4C1C2+1)C3 
^5 ^ 4C2(C3-Ci)-1- 



(137) 



2. NOR with live lapse and shift 



The most general conserved energy can be written as 

10 
£ = ^c,e„ (138) 



with eo, ..., £5 as given before. Defining 



sX, + {bs-^)U + il-bs)^j 



£S 



Y. + ^ik) 



Zi = fis 



sX,-H(fe-^)| 



£5 



Y, + ^U 



(140) 



we have five new energy terms containing derivatives of 
the shift vector 



£6 = P'J'.,, 

£7 = P\^P\, - i3\,i3\,, 

£9 = K'^p\,-'^W'd,, 



£10 = KP' 



-zx. 



(141) 
(142) 
(143) 
(144) 

(145) 



Assuming 6 7^ and e^ 7^ ^l (two necessary condi- 
tions for strong hyperbolicity in NOR), the coefficients 
C3, C4, C5 are always determined by the rest of the coeffi- 
cients: 

4&C3 = 2aco -I- 4s/isC8 -I- (2a-f s/is)c9, (146) 
(£l-Ml)c4 = 2(a + d)co-H2(a-|-3d)ci-|-2s/isC6(147) 
+ {a + d + sfis)cg + {a + 3d+ s/is)cio, 

65-1 
(£i-ML)c5 = £5C6 - (co + Ci) H (Cg -h CioX-148) 
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Energy conservation imposes 4 more conditions: 

{ab-fis)cg + 'ifis{bs-l)cs = 0,(149) 
2eLCQ + [2ab -l-2c + 2bd- elcs 
+pLs{p - 2)]cg + 2{ab - 1 - 3c + 36d)cio 

+4/xs(6s - 1)C6 = 0,(150) 

[2 - 2ab + &c - Qbd - ELes 

+^J■L + usip- 2)]cio 

+2[fis{p - 2bs) - esHL]ce 

+2eLCi + {l-2ab+2c-2bd + ^iL)c9 = 0,(151) 

co(l - 2ab + ^i + 2c - 2bd) 

+2ci(l — a6 + 3c — 3bd) + /^^(l — bs)cio 

+ l[^s{p-2bs)-esfiL]c9 - 0.(152) 

For generic values of the parameters, these four condi- 
tions allow us to solve for cg, cs, cg, cio, leaving only four 
free coefficients co,ci,C2,C7. It is possible, however, to 
choose the parameters so that some of those four condi- 
tions are automatically obeyed, and then the number of 
free coefficients is larger. For example, for NOR with im- 
plicit harmonic gauge and s = 1 all four conditions are 
automatically obeyed, leaving eight free coefhcients c^. 
Note that the condition (|134|l has now been transformed 
into H152|l . Note also that for e^ = conditions (|149l 
I150lll5l|l only involve the new coefficients ce, ..., ciq. 

So far we have only looked at energy conservation. 
Symmetric hyperbolicity requires positivity of the con- 
served energy. Ideally, we would like to translate pos- 
itivity into two sets of inequalities: ranges of evolution 
parameters that lead to a symmetric hyperbolic system, 
followed by (parameter-dependent) inequalities for the 
free coefficients of the conserved energy given above. We 
achieved this for NOR with fixed gauge in i^], but have 
not managed for the systems studied here. 

Sometimes, for example for the LS system 8], it is pos- 
sible to choose the coefficients of a generic energy to make 
it positive, and then solve the energy conservation equa- 
tions for the evolution parameters of the system. This 
method proves that the system is symmetric hyperbolic 
for some values of the parameters, but works only if there 
are enough parameters in the evolution system. It may 
also be difficult to find energies relating to specific inter- 
esting values of the evolution parameters. 

For NOR with LS gauge, we only have 12 parame- 
ters, namely (a, 6, c, d), {fiL,Ps, £l, £s) and (p, p, a, s), to 
solve 15 linear equations imposed by energy conserva- 
tion on 19 coefficients of a general energy. Of the 19 
coefficients, 11 multiply squares (diagonal terms) and 8 
multiply mixed products (off-diagonal terms). Setting 
7 of the off-diagonal terms equal to zero, we have been 
able to solve the remaining positivity conditions and all 
energy conservation equations to obtain a 9-parameter 
(7 evolution parameters and 2 energy coefficients) fam- 
ily of symmetric hyperbolic systems. Three of the five 
relations among the parameters are e^ = 0, 6s = 1 and 



Pl^S = (p — 2)/zs; the other two are rather longer. How- 
ever, there is no reason to assume that any of the positiv- 
ity conditions fail if all off-diagonal coefficients of the en- 
ergy are given sufficiently small non-zero values. There- 
fore we believe that our 9-parametcr family is embedded 
in a 16-parameter open set of symmetric hyperbolic for- 
mulations, which could be expressed in terms of the 12 
evolution parameters and 4 energy coefficients. 

In particular, we have been able to construct positive 
and conserved energies for relevant choices of the param- 
eters: for example for NOR-B with /x/, = /i5 = s = 1, 
El = 0, £5 = — 1, the four conservation equations reduce 
to cg{p— 1) = and (2c6-l-cio)(p— 1) = 0. For p=\ the 
eight coefficients are free and a possible positive energy 
is given by 



Co = 1, ci = 112, 


C2 = 


3 

16' 


C6 = 12, 


(153) 


C7 = — 5 cg = 16, 


cg = 


-8, 


Cio = -8. 


(154) 



We can also construct a positive energy for BSSN-C (p — 
p = 2/3, a — 0) choosing 6 free coefficients, although we 
have only been to solve the inequalities by finding specific 
sets of energy coefficients, rather than by giving ranges. 



V. HYPERBOLICITY WITH IMPLICIT 

HYPERBOLIC LAPSE AND EXPLICIT 

HYPERBOLIC SHIFT DRIVERS 

Recently, Campanelli et al. [13 and Baker et al. [20J . 
followed by Diener et al. |2l| et al. and Herrmann et 
al. [231 have reported significant progress in binary black 
hole evolutions with the BSSN-C formulation using an 
implicit hyperbolic lapse driver combined with an explicit 
hyperbolic shift driver. The gauge choice differs slightly 
between these four groups. We shall look at a family that 
includes all four cases plus a fifth that we suggest as a 
simpler alternative. In the following we look at BSSN- 
C in the form of NOR-B {a = b=l, c = d= -1/3) 
with the further choice s = I, p = p ^ 2/3, so that 
j« — 7i/'^r*. We also set ct = 0. The gauge choice we 
consider is 



{do + CLd^)lna c 


- -PlK, 


(155) 


{do + Csd^)P' c 


- b\ 


(156) 


{do + Cbd^)b' z 


- psa{do + Cfd^) f. 


(157) 



(In each case, the factor of proportionality between our 
variable 5' and the variable B^ in which these gauges are 
originally formulated is different. We have introduced V 
to simplify the comparison of the principal terms.) In 
the principal part, this corresponds to the gauge choices 
of the four groups with the values of the six parameters 
given in Table Q] All four groups set (s — (1, = 1. By 
eliminating the variable 6* in favour of a second time 
derivative of /3*, and noting that 6* does not appear 
anywhere else in the evolution equations, we see that 
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iPsjPb) = (0:1) and (1,0) would be equivalent in the 
principal part. 

With the shorthands 



(158) 

(159) 
(160) 
(161) 



the second column of Table Q] then gives all characteristic 
speeds of these four systems. In all four systems strong 
hyperbolicity breaks down when some of the speeds co- 
incide. The conditions for avoiding this are listed in the 
third column of the table. It is possible that these in- 
equalities hold for the specific binary black hole evolu- 
tions because /3' was sufficiently small in them, but one 
would expect numerical difficulties to appear for /3' suf- 
ficiently large to violate one of these conditions, for ex- 
ample in co-rotating coordinates. It would be simple and 
interesting to verify this. 

If we set (i^ — (g = (i, = Q = 0, all speeds become 
independent of /3* . They are 



In the resulting linear problem with constant coeffi- 
cients we take a Fourier transform 




(0, ±1, ±^jiE, ±^jr^, ±v/Vs/3) 



(162) 



and the conditions for strong hyperbolicity are /ig ^ /i^ 
and ^s 7^ 3/Xi/4, which can now easily be arranged 
through the choice of /j-l and ^s a-s functions of a, for 
example /^s — 3/4 and ^l = 2/q!, which would not coin- 
cide as long as a < 2, independently of the shift. Strong 
hy perb olicity of BSSN with this gauge choice was shown 
in [23 ■ Hyperbolicity of the Z4 formulation with a fam- 
ily of gauges with the same principal part and, in our 
notation, C/ = Cs = and Cl and C,b free, was examined 
in il. 



VI. MODE ANALYSIS RESULTS 

A. Formalism 

If the complete evolution system is not hyperbolic, we 
can use mode analysis to obtain at least necessary condi- 
tions for well-posedness. We make two approximations: 
we linearise around a background, and we approximate 
all background-dependent coefficients of the linearised 
equations as constant ("frozen"). Physically this corre- 
sponds to investigating small amplitude, high frequency 
perturbations. (The results obtained will in general de- 
pend both on the background solution, and the choice of 
gauge on it.) 



u{x\t) 



u{uJi, t) d LO 



(163) 



and consider the evolution of one Fourier mode at a 
time. Well-posedness of the Cauchy problem means that 
||'"(-,i)|| < /(i)||'u(-,0)|| with/(i) independent of the ini- 
tial data u(a:*, 0). This is the case if and only if a similar 
bound can be obtained for u{uJi, t) with f{t) independent 
of u)i. To establish well-posedness it is sufficient to con- 
sider only certain leading order terms. This is well-known 
for strongly hyperbolic first-order systems 26], but we 
need a generalisation which is given in the following 

Theorem 1: Consider a linear problem with constant 
coefficients which in pseudo-differential form is 

dtu{uj„t) = [M{Lj,) + M'{uj,)]u{LU,,t), (164) 
M(wi,0) = uo{uji). (165) 

The initial-value problem is well-posed if M{uJi) is diag- 
onalisable, 



M{ujf)=T-\uj,)K{uj,)T{uj,) 



(166) 



with A(aji) diagonal, and if M{uJi) and M'{uJi) obey the 
bounds 



\T{oj,)\\T-\uji)\ < K,, 
|ReA(c^,)l < K2, 



(167) 
(168) 
(169) 



with Ki , K2 and K^ independent of tOi . 

For a more general statement in the context of semi- 
group theory see [25J . An elementary proof of our version 
of the theorem is given in Appendix ^ and a technical 
complication in applying the theorem to second-order in 
space systems is resolved in Appendix 151 

Clearly the well-posedness of the linearised, frozen- 
coefficients problem is necessary for the well-posedness 
of the full problem. For strongly hyperbolic problems 
well-posedness of the linearised frozen-coefficient prob- 
lem around any background, together with smoothness 
of T{uJi), is also sufficient _26j, but we do not know if 
that is true for the wider class of problems considered 
here. 



B. Implicit hyperbolic drivers 

To establish notation we recast the problem of strong 
hyperbolicity of the NOR formulation with implicit hy- 
perbolic lapse and shift drivers, which was already dis- 
cussed in Section IIVI into the language of pseudo- 
differential operators. We define 



■u = {iujv,w), 
Z = ibjz. 



(170) 
(171) 
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Cl Cs Ci> C/ ^^L 



fiS 



speeds 



strongly hyperbolic for 



Campanelli et al. 
Baker et al. 
Diener et al. 
Herrmann et al. 
Beyer-Sarbach 



1 1 1 
1 1 



|a-^y/^ 



1111 ^^b"l 3 gP-^ 1/3 
-^ -^ -^ -^ a 4 VP- ' 



1110^ 

a 

/ 



fBL 

.-271/3 

GH 



(0, ±1, i^/ir, -/3„, A2±, ±^/7Ii^ 

(0, ±1, ±^/JlL, Asi, A4±) 

(0,±l,-^„,Ai±,A2±,±yM?) 

(0, ±1, — /3n, Ai±, Aai, A4±) 

(0, ±1, zfcVMT, zfc^M?, ± yi^WS) 



1^1 < K 3ml- 4 ^5)7(3^^1) I 

Mi / 4/is/3, i/3! < ^A^ 
|^|<|(3ML-4Ms)/(2y37^) 
MS / MJ-) MS / 3mi,/4 



\/4ms/3| 



TABLE I: Parameter values, characteristic speeds, and conditions for strong hyperbolicity for the four "puncture evolution" 
codes and the variant of Beyer and Sarbach. For comparison with these authors, we have expressed our parameters ^l and 
/is in their notation, but for comparison between the gauges we give the characteristic speeds in our notation. 



where w = a/t^^w^Wj and z = (lna,/3*). The factors 
ibj appear because this is a pseudo-differential reduction 
to first order of a second-order in space system Q- We 
define the pseudo-differential equivalent of the derivative 
operator 9o, 



use a different pseudo-differential reduction of the form 



dn=a ^{dt~ ioj(3nI). 



(172) 



The evolution of the gauge variables has the principal 
part z ^w + dv + dz, and combining this with (|2I3(I the 
coupled pseudo-differential equations are 




P Q 

V W 




(173) 



where P, Q, V and W are independent of lo but depend 
on Hi = uji/oj. Here ~ means that terms corresponding to 
the lower-order part M'{ijji) have been neglected. (In the 
context of a pseudo-differential reduction to first order 
there is a slight complication in the definition of M' [uji)^ 
which is discussed in Appendix IBI) P and Q can be read 



off from (|2()I22(I . The nontrivial part of the matrix V is 
given by l|S^ . W contains the parameters Cl, Cs? ^l and 
es defined in (|111I112|I . 

We have analysed this system above in Sect. IIVI and 
do not discuss it further here. 



C. Explicit parabolic drivers 

These have the general form z ^ w -\- dv + ddz. The 
pseudo-differential form of the coupled equations is 



iuP iujQ 

'Uj-R + iujR' —uj-S + iujS' 





(174) 



where the matrices i?, S*, R' and S" are independent of 
LO. Note that in order to apply Thm. 1 we need to work 
with the lower-order terms R' and S' . To calculate them 
we need to explicitly linearise the gauge conditions. If 
the principal part matrix defined in 1174(1 obeys the con- 
ditions of Thm. 1, we obtain a necessary condition for 
well-posedness in L^{w,dv,dz). Alternatively, we can 




iujP —uj'^Q + iujQ' 
iujR ^Lo'^S + iLoS' 




(175) 



The conditions for Thm. 1 are then a necessary condition 
for well-posedness in Lp'{w, dv, z). The lower-order terms 
that must be kept are now S" and Q' . 

When either set of lower-order terms is kept, the re- 
sulting principal part matrix M{u!i) becomes very com- 
plicated. In particular, the eigenvalues depend on the 
background spacetime, and the matrix no longer splits 
into block-diagonal form with scalar, vector and tensor 
blocks. In order to make progress, we linearise around 
Minkowski spacetime using Cartesian coordinates. 

In the Cartesian Minkowski case the lower order terms 
i?', S' and Q' all vanish, and we are left with the simpler 
matrices M{uji) of the form 



itoP iujQ 




(176) 



which are similar, and therefore equivalent for the pur- 
poses of Thm. 1. The diagonalisability and eigenvalues 
of these matrices already give interesting necessary con- 
ditions for well-posedness. 

We consider explicit parabolic drivers in the form 



9o In a 






(177) 
(178) 



Note that when only principal terms arc considered we 
can pull the factor of a out of the time derivative in 
(|178|l . R and S can then be read off from ||Hlll and ^^. 
We find that the simplified necessary conditions based on 
linearisation around Minkowski are obeyed for this gauge 
with ADM and NOR-B, but not NOR- A (the scalar sec- 
tor fails to be diagonalisable) , for p < 2, i^l > 0, i^s > 
and t^i ^ (2 — p)vs- The eigenvalues are 0, zLiuj, — i/^o;^, 
and —(2 — p)vs^'^- 

The lapse condition H177|l with fixed shift together with 
ADM or NOR-A/B is ill-posed according to Thm. 1. This 
can be understood roughly as follows: The principal part 
of the lapse evolution equation is d ~ Aa, which is au- 
tonomous, a can therefore be considered as a a given 
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function in the principal part of the remaining evohition 
equations, and so their principal part is the same as for 
fixed (undensitised) lapse (and fixed shift). 

To investigate the explicit parabolic drivers using dt 
on both the left and right-hand sides, that is, (|92I96() . we 
consider the matrix 



iuP 



u^R + iuPnT -uj'^S 




(179) 



(Again we drop the lower order terms by linearising 
around Minkowski.) Neither ADM or NOR-A/B_are di- 
agonalisable with this gauge, except of course for /3„ = 0. 
Because the explicit parabolic lapse driver with fixed 
shift using do alone is not well-posed (as we saw above), 
neither is the version using dt as the two coincide for 

(3n = 0. 



D. Implicit parabolic drivers 



These have the general form z ~ 
sation, in pseudo-differential form 



w^ 



{iLuT + T')u, 



-dv. After lineari- 



(180) 



where T and T' are independent of w. In order to calcu- 
late the lower order term T', we again need to explicitly 
linearise the gauge conditions. Eliminating Z, we find 



dou' 



-uj^QT + iLjiP + QT')] 



(181) 



We can again make progress by linearising around 
Minkowski spacetime in standard coordinates, so that 
T' vanishes. T is equal to V with /i^ and /ig replaced by 
vl and ^^5. 

The conditions of Thm. 1 are obeyed for the implicit 
parabohc drivers (ESEZI) with ADM and NOR, with the 
same inequality conditions as for the explicit parabolic 
drivers. The eigenvalues are also the same as for the 
explicit parabolic drivers. 

The implicit parabolic lapse driver H93|l with fixed shift 
obeys the conditions of Thm. 1 for NOR-B for i/^ > 0. 
The eigenvalues are 0, ztiw and —i'luj'^. This is not true 
for NOR-A (because the scalar sector is not diagonalis- 
able) or ADM (because the vector sector is not diagonal- 
isable) . 



E. Elliptic freezing conditions 



These have the general form w + dv : 
sation, in pseudo-differential form. 



0. After lineari- 



~ {-uj^R + iujR')u + {-Lu'^S + iLuS')Z, (182) 

where R', R, S' and S can be chosen to be the same 
matrices as above. Eliminating Z, we obtain the reduced 
system 



dou ~ iuj{P - QS~'^R)u. 



Note that this algebraic reduction can only be carried 
out in pseudo-differential form. Intuitively speaking, R' 
and S" do not appear in the principal part because of the 
smoothing effect of solving an elliptic equation. 

The elliptic freezing conditions doK = 0, dof^ — 
obey the conditions of Thm. 1 together with ADM and 
NOR-B but not NOR-A (the scalar sector is not diago- 
nalisable) for p ^ 2. The speeds are and ±1. 

The elliptic lapse condition alone, with fixed shift, is 
not diagonalisable for either of the three systems. This 
slightly surprising result can be explained as follows: The 
only component of doKij ~ — (In a) _ jj + . . . in the pseudo- 
differential approach is doKnn — a;^(lnQ;) -I- . . .. But the 
elliptic equation for a is — u;^(lna) ~ 0. The scalar sec- 
tor of the principal part is then the same as for fixed 
(undensitised) lapse, and so is not diagonalisable. 

However, if one formally enforces the constraint K = 



(x) by setting Kq 



-Knn the scalar sector of NOR- 



A/B (but not ADM) becomes diagonalisable, and the 
resulting system obeys the conditions of Thm. 1. In the 
BSSN formulation K is an explicit variable, and one can 
enforce K = (^{x) by not evolving this variable when 
BSSN is used with maximal slicing and fixed shift |27l| . 
With this modification the system is well-posed, but if K 
was evolved it would be ill-posed. 



Implicit hyperbolic drivers with parabolic 
damping terms 



The gauge conditions H1U5I1U8|I . linearised around 
Minkowski in Cartesian coordinates, give 




iojP iLoQ 

iLuV~2uj^R iujW-2Lu^S 




(184) 



For simplicity we restrict the hyperbolic gauge conditions 
to ( = €l = es ^ 0, and we consider ADM. Without the 
damping terms, that is with vl = ^s = 0, the system 
is strongly hyperbolic for /i^ > 0, /is > and p < 2, 
with eigenvalues 0, ±zw, ±iujy/]IZ, ±iu}^{2 — p)ps and 
zLiLOy/JTs, provided that fj,L 7^ 1, ps 7^ 1/(2 — p) and 
Ps ¥" Ml/(2 — p), corresponding to real speeds. With 
the damping turned on, the pair of eigenvalues ±iuj^/JJLL 
becomes 



UL ± \luj^v'i 



^^PL, 



(185) 



and similarly for the other eigenvalues. This has neg- 
ative real part for all oj > 0. However, the scalar sec- 
tor fails to be diagonalisable if these two eigenvalues 
coincide, that is for to — ^ plI^l- Similar problems 



(183) 



arise for u: — y/pZ/{^/2 — pvl) in the scalar sector and 
uj ::= J PlI^l in the vector sector. In numerical appli- 
cations, choosing vg^vj^ oc Ax with a sufficiently small 
constant of proportionality makes the problematic wave 
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numbers larger than the grid frequency, and this ill- 
posedness should then not lead to numerical problems 
\ll\ . Similar results hold for NOR-A/B at least in the 
special cases s = or p = 2. 



VII. CONCLUSIONS 

Continuing our research programme on second-order 
in space, first-order in time formulations of the Einstein 
equations y, S l^^i we have reviewed the relationship 
between the BSSN formulation of the Einstein equations 
that is widely used in numerical relativity, and the NOR 
formulation (NOR-A) that was suggested as a simpler 
alternative to BSSN. We have showed that a variant of 
NOR, (NOR-B) has the same principal part as BSSN 
when the latter is formally restricted to solutions of the 
algebraic constraints (BSSN-C). 

We have shown that the principal parts of the KST for- 
mulation and the NOR formulation are equivalent when 
the gauge is fixed. Each formulation has five parameters, 
and we have given four relations between the two sets 
of parameters. One parameter in each formulation does 
not have a counterpart in the other formulation and does 
not influence the level of hyperbolicity. The analogy also 
holds with an evolved lapse, and symmetric hyperbolic 
cases exist for both fixed gauge and live lapse. With an 
evolved lapse and shift, KST has four parameters more 
than NOR with the same gauge. 

We have reviewed various differential equations for the 
lapse and shift currently in use as coordinate conditions 
in numerical relativity, with an emphasis of their origin 
as symmetry-seeking coordinates, and their classification 
as elliptic, parabolic, hyperbolic, or neither. These are 
all based on X-freezing shift and F-freezing lapse, which 
in turn are related to the well-known maximal slicing, 
minimal distortion gauge. 

From the point of view of "symmetry-seeking" gauge 
conditions it is more natural to have evolution equa- 
tions for the gauge of the form dta,dt/3 — . . ., while 
the ADM evolution equations are naturally of the form 
i9o7y", d^Kij — . . .. With this mixture, there are two sets 
of light cones in the system, one centred around n° and 
one around id/dt)"'. Where these overlap, there is a dan- 
ger of strong hyperbolicity breaking down, but we have 
shown that some dt gauges can be implemented in do 
form, which would avoid this problem. 

We have analysed the hyperbolicity of the ADM, NOR- 
A and NOR-B formulations with the most general im- 
plicit hyperbolic lapse and shift drivers. Interestingly, 
NOR-A and NOR-B can be made strongly hyperbolic 
even when using some live gauges with dt- 

We have also investigated the hyperbolicity of NOR-B 
(equivalent to BSSN-C) with a family of implicit hyper- 
bolic lapse and explicit hyperbolic shift drivers that in- 
cludes the gauges used by four different research groups 
in "moving puncture" evolutions. We find that all these 
gauges become ill-posed for large enough values of the 



type 


lapse -1- shift 


lapse only 


impl. hyp. using do 


ADM, NOR-A/B 


NOR-A/B 


impl. hyp. using dt 


NOR-A/B 


NOR-B 


expl. par. using do 


ADM, NOR-B 


no 


expl. par. using dt 


no 


no 


implicit parabolic 


ADM, NOR-A/B 


NOR-B 


elliptic 


ADM, NOR-B 


no 


impl. hyp. using do 


ADM, NOR-A/B 


? 


with par. damping 







TABLE II: Overview of well-posedness results for combina- 
tions of gauge conditions and formulations of the Einstein 
equations. In the upper half of the table, strong hyperbolic- 
ity has been checked. In the lower half of the table only only 
a necessary condition (Thm. 1) has been checked. A question 
mark indicates that we have not done the calculation. 



shift, and propose a simple modification (replacing dt by 
do) that is strongly hyperbolic for arbitrary shift. 

For certain families of elliptic gauge conditions, 
parabolic drivers, and hyperbolic drivers with a heat 
equation-type damping, where the notion of hyperbol- 
icity does not apply, we have carried out the mode anal- 
ysis and have checked a necessary condition for well- 
posedness of the Cauchy problem. 

Khoklov and Novikov 28] have investigated the well- 
posedness of a number of gauges independently of any 
particular formulation in the frozen coefficient approx- 
imation. Their well-posedness is therefore a necessary 
condition for ours. We review their method and compare 
results in Appendix lUl 

Our results for puncture gauges are summarised in Ta- 
blelHand our results for all other gauges in Table UTI They 
are positive for most combinations of gauge and formu- 
lation, with a few interesting exceptions. They provide 
a theoretical underpinning to some conditions already in 
use in numerical relativity, and suggest certain improve- 
ments. 

All equations in section IIVI were derived using 
xTensor, an open-source Mathematica package for ab- 
stract tensor calculations, developed by J.M.M. It 
is available under the GNU Public Lice nse from 
http://metric.iem.csic.es/Martin-Garcia/xAct 
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APPENDIX A: PROOF OF THEOREM 1 

We go to the diagonal basis U — Tu. It obeys 

dtU = (A -I- TM'T-^)U, (Al) 
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and so its norm squared obeys 

dtiU^U) < 2(^2 + K,K3)U^U. (A2) 

This can be integrated to give 

|?7(a;.,i)|<e(^^+^i^^)*|f7(o.„0)|, (A3) 

and going back to the original basis we have 

\ui^.,t)\ < Xie(^^+^^^^)*|zi(c.„0)|, (A4) 

that is, the growth is bounded independently of the 
wavenumber cu. Taking the inverse Fourier transform and 
using Parseval's Theorem, we obtain the L^ estimate 

||w(-,t)||<i^ie(^=+^^^^)*||w(-,0)|| (A5) 

in real space. 

The well-known result that lower-order terms do not 
affect the well-posedness of strongly hyperbolic systems 
(see for example Theorem 4.3.2 of [2^) is a special case 
of Thm. 1. A system of linear first-order evolution equa- 
tions is strongly hyperbolic if M{uji) = iojiP"^ = iujPn, 
with uji = Lurii and P„ = UiP'', and where P„ is di- 
agonalisable with real eigenvalues for every Ui (so that 
ReA = 0), T depends smoothly on n^ but not on w (and 
so is bounded) , and M' is independent of tOi (and so is 
bounded). 



APPENDIX B: PSEUDO-DIFFERENTIAL 
REDUCTION TO FIRST ORDER 

The pseudo-differential form of the linearisation of 
(I2I3|) on a constant background, neglecting the terms in 
z and writing all (linearised) lower-order terms, is 




iLuA + A' 
-u?C + iujC 



C" iujD 




where the matrices A, A' etc. are independent of to but 
depend on rii . Replacing w by 1/^ = itov multiplies the first 
row by iuo and divides the first column by iu, so that the 
highest power of w is the first power. This constitutes a 
pseudo-differential reduction to first order. 

The terms proportional to iuj form M{uJi) and the 
terms of 0{\) go into M'{ijji). The term (icj)^^C" poses 
a problem, as it is not bounded independently of uJi . In- 
tuitively, it should be part of M'{LOi) because we are in- 
terested mainly in the limit u; — > cx). We can obtain a 
rigorous result to the same effect if we retain v in the 
system. We then have 





For the purposes of Thm. 1 we only need to analyse 
M{ijji), which here only means the block {{A, B), (C, D)), 
which is called P in the main text. 



APPENDIX C: ANALYSIS OF PURE GAUGE 
PERTURBATIONS 



The effect of a set of gauge conditions can be evaluated 
in purely geometric terms, independently of any formu- 
lation of the Einstein equations. In this approach, one 
considers a fixed spacetime obeying all 10 Einstein equa- 
tions. One evolves only the four coordinates x^^ and the 
vector {d/dt)°- on this spacetime background, treating 
x'^ as four scalar fields and {d/dty as a vector field on 
the background spacetime. To make connection with the 
3-1-1 split, let n"" be the unit normal on the surfaces of 
constant a:°, and let {d/dt)°- = arf'+fi"-, where /3°na = 0. 
The resulting non- linear evolution system for (x^ .a. ff^) 
and a numerical implementation are described in |30j. 

A linearisation of this scheme is used in |23| to inves- 
tigate the well-posedness of gauges. The gauge pertur- 
bation is parameterised by a vector field V'a, which gives 
rise to a perturbation of the 4-metric 5gab ~ ^"^^ la^b)- 
The 10 components of this equation can be written in a 
3+1 split as 



(CI) 

(C2) 
(C3) 



(B2) 



^0 = a5a-P'5l3,-l3'(i'i,,,, 
(57„- = -2V'(,j) + 2r^„?A^- 



5Kij can be obtained by taking a time derivative of (|C3p 
and using the definition of the extrinsic curvature tensor 
(O. These equations are coupled with algebraic, elliptic 
or evolution equations for 5a and 5j3i in terms of 5^ij and 
5Kij . Clearly the well-posedness of any gauge defined in 
this way is a necessary condition for that gauge to be 
well-posed together with any specific formulation of the 
Einstein equations. 

28] shows that the combination a = eK, /?* = (ex- 
plicit parabolic ii'-driver with zero shift) is ill-posed. By 
taking the limit e = 0, it is also claimed ivT-freezing lapse 
with zero shift is ill-posed, although this limit is singular. 
Our results agree with these claims. However, we only 
consider the linearisation around Minkowski spacetime 
in Cartesian coordinates, which simplifies the calculation, 
while in 28J all lower-order terms in an explicit linearisa- 
tion are kept, and are found to give an unbounded ReA 
for the parabolic lapse driver with fixed shift. On the 
other hand we find that M{uji) is not diagonalisable when 
perturbing around the Minkowski background, while di- 
agonalisability is not checked in p8| . 
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